QCD quark condensate in external magnetic fields 
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We present a comprehensive analysis of the light condensates in QCD with 1+1+1 sea quark flavors 
(with mass-degenerate light quarks of different electric charges) at zero and nonzero temperatures of up 
to 190 MeV and external magnetic fields B < 1 GeV 2 /e. We employ stout smeared staggered fermions 
with physical quark masses and extrapolate the results to the continuum limit. At low temperatures we 
confirm the magnetic catalysis scenario predicted by many model calculations while around the crossover 
the condensate develops a complex dependence on the external magnetic field, resulting in a decrease of 
the transition temperature. 
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Introduction. — Strong (electro)magnetic fields promi- 
nently feature in various physical systems. They play an es- 
sential role in cosmology, where magnetic fields of 10 14 T 
and 10 19 T may have been present [1,2] during the strong 
and electroweak phase transitions of the universe, respec- 
tively. Magnetic fields with strengths up to B ~ lfj 14-16 T 
(V eB ~ 0.1 — 1.0 GeV) are also generated in non-central 
heavy ion collisions [3-6] at the Relativistic Heavy Ion 
Collider (RHIC) or the Large Hadron Collider (LHC). Fur- 
thermore, for certain classes of neutron stars like magne- 
tars, magnetic fields of the order of 10 10 T have been de- 
duced [7]. In addition to this phenomenological relevance, 
external (electro)magnetic fields can be used to probe the 
dynamics of strongly interacting matter, i.e. the vacuum 
structure of Quantum Chromodynamics (QCD). 

One of the most important aspects of QCD is chiral sym- 
metry breaking. At zero quark masses the chiral conden- 
sate ipip is an order parameter. It vanishes at high tem- 
peratures where chiral symmetry is restored but develops a 
nonzero expectation value in the hadronic phase. In nature, 
quark masses are nonzero and the corresponding quark 
condensates, though only approximate order parameters, 
still exhibit this characteristic behavior around the transi- 
tion temperature between the hadronic and the quark-gluon 
plasma phases. Lattice simulations revealed that for physi- 
cal quark masses this transition is an analytic crossover [8], 
leading to a transition temperature T c which may depend 
on the observable used for its definition. 

In the response of QCD to external magnetic fields, 
'magnetic catalysis' refers to an increase of the conden- 
sate with B. This implies a S-dependence of T c as 
well. Almost all low-energy models and approximations 
to QCD [9-35] as well as lattice simulations in quenched 
theories [36, 37] and at larger than physical pion masses in 
N f = 2 QCD [38, 39] and in the N f = 4 SU(2) theory 
[40] found ipip(B) and T C (B) to increase with B. Ex- 
ceptions in this respect with a decreasing T C (B) function 
are the results obtained within two-flavor chiral perturba- 



tion theory [41], in the linear sigma model without vacuum 
corrections [42] and in the bag model [43]. 

In contrast to the majority of the above results, our large- 
scale study of QCD in external magnetic fields with phys- 
ical pion mass M„ = 135 MeV and results extrapolated 
to the continuum limit [44] has revealed the transition tem- 
perature to decrease as a function of the external magnetic 
field. This applies to the T c 's defined from the quark con- 
densate, the strange quark number susceptibility and the 
chiral susceptibility. In particular, we found the condensate 
to depend on B in a non-monotonous way in the crossover 
region. 

In Ref. [44] we have already pointed out two rationales 
why former lattice studies are at variance with these recent 
findings: coarser lattices and larger quark masses. Obvi- 
ously, it is also very important to address the differences 
between our QCD results and many model and chiral per- 
turbation theory (xPT) predictions, especially since the lat- 
ter methods can be used to investigate regions that are not 
easily accessible to lattice simulations, e.g., QCD at a non- 
vanishing baryon density. 

In this Letter, we present a detailed analysis of the de- 
pendence of the light quark condensates on B and on the 
temperature T, based on the T > simulations described 
in [44] and new simulations at T = 0. The data — all 
continuum extrapolated — are presented in ways that will 
enable to refine model assumptions and parameters. We 
also perform a first comparison to x?T and to Polyakov- 
Nambu-Jona-Lasinio (PNJL) model predictions. We aim at 
a better understanding of the physical mechanisms behind 
the differences. This in turn should be of phenomenologi- 
cal relevance. 

Below we introduce our notations and the simulation 
setup. We then present our results and compare them to 
XPT and PNJL predictions, both at zero and nonzero tem- 
peratures. 

Condensate on the lattice at nonzero B. — We study 
QCD coupled to a constant external magnetic field B, 
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pointing in the positive z direction. Such a field can be 
implemented by multiplying the U £ SU(3) links of the 
lattice by complex phases. The specific choice of these 
phases and our setup are detailed in Ref. [44]. The exter- 
nal field couples only to the quark electric charges qf with 
/ labelling the different flavors. Thus, the magnetic field 
only appears in combinations qfB. 

In a finite periodic volume, the magnetic flux is quan- 
tized [45, 46]. This quantization on a lattice with spacing 
a amounts to, 



(N s a) 2 • q d B = 2nN b , N b e Z, 



(1) 



where the smallest quark charge (that of the down quark), 
\Qd\ = e/3 enters, with e > being the elementary charge. 
Here N s is the number of lattice sites in a spatial direction 
(our lattices are symmetric in space). Similarly, N t counts 
the lattice points in the temporal direction. The spatial vol- 
ume of the system is given by V = (N s a) 3 and the temper- 
ature is related to the inverse temporal extent of the lattice 
asT= (N t a)-\ 

The quark condensate can be derived from the partition 
function, which in the staggered formulation of QCD with 
three flavors (/ = u, d, s) is given by the functional inte- 
gral, 
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where f3 = 6/g 2 is the inverse gauge coupling, S g the 
gauge action and M(U, qB, m) = lp(U, qB) + ml the 
fermion matrix. For S g we use the tree-level improved 
Symanzik action, while in the fermionic sector we employ 
a stout smeared staggered Dirac operator ij). The details of 
the lattice action can be found in Refs. [44, 47]. The lattice 
sizes range from 24 3 x 32 to 40 3 x 48 for the zero temper- 
ature simulations, while at non-vanishing T we investigate 
24 3 x 6, 24 3 x 8 and 28 3 x 10 lattices. We set the quark 
masses to their physical values, with mass-degenerate light 
quarks: m u = m d = m ud . The electric charges of the 
quarks are q d = q s = —q u /2 = —e/3, therefore we 
need to treat each flavor separately. The line of constant 
physics (LCP) [m ud (/3),m s (/3)] was determined by fix- 
ing the ratios M^/f K and M K /f K to the experimental 
values. The lattice spacing a(f3) is defined by keeping 
fx = f% (P)/ a (P) fixed, for details see Ref. [48]. At 
T = the continuum limit a — > corresponds to (3 — > oo. 
At nonzero temperature, it is convenient to define the con- 
tinuum limit as N t —> oo, keeping T fixed. 

The quark condensate is defined as the derivative of In Z 
with respect to the lattice mass parameter 

- , ,„ m . T d\nZ(B, T) 
W> f (B,T)==- Q V (3) 
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To carry out the continuum limit, the lattice condensate 
•ipil) needs to be renormalized since it contains additive (for 
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FIG. 1, The change of the renormalized condensate due to the 
magnetic field at T — as measured on five lattice spacings and 
the continuum limit. 



m > 0) and multiplicative divergences. These cancel [44] 
in the following combination, 

2m„ d 



J u,d 



[B,T) 



M1F 2 



[i;^ d (B,T) - ^ d (0,0)]+l, 

(4) 

where, to obtain a dimensionless quantity, we divided by 
the combination M 2 F 2 which contains the zero-field pion 
mass M n = 135 MeV and (the chiral limit of the) pion de- 
cay constant F = 86 MeV [49]. This specific combination 
enters the Gell-Mann-Oakes-Renner relation, 



2m ud -W(0,0) = M 2 F 2 + 



(5) 



Note that the normalization in definition (4) can easily be 
converted into the slightly different one employed in for- 
mer studies by the Budapest- Wuppertal collaboration (e.g. 
Refs. [8, 50]) and in Ref. [44]. We define the change of the 
condensate due to the magnetic field as 

AZ u . d (B, T) = E U)d (S ; T) - E Uid (0, T). (6) 

Note that the ^ip(0, 0) term cancels from this difference. 
In our normalization, Eq. (6) defines the change of the con- 
densate caused by a nonzero B, in units of the chiral con- 
densate at B = and T = 0. This normalization will be 
advantageous when comparing the lattice results to xPT 
and model predictions, which are usually given in units of 
ipip(0, 0). The +1 is included in Eq. (4) so that the chiral 
limit of the condensate is fixed to 1 at T = B = 0, and ap- 
proaches as T — > oo. At nonzero quark mass S U)( j will 
still start from 1 at T = B = 0. At very high tempera- 
tures, however, it is well known from the free case [51, 52] 
that the condensate receives a contribution ~ mT 2 . This 
term is negligible for the temperatures under study and it 
cancels exactly from AS U d . 

Results. — In Fig. 1 we display the renormalized differ- 
ence A(S U + Ti d )/2 as a function of B at T = 0, for five 
different lattice spacings. We carry out the continuum limit 
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FIG. 2. Continuum extrapolated lattice results for the change of 
the condensate as a function of B, at six different temperatures. 
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FIG. 3. Comparison of the continuum limit of the change of the 
condensate to the xPT [13-15, 54] and the (P)NJL model [18, 55] 
predictions. 



by fitting the results to a lattice spacing-dependent spline 
function (for a similar fit in two dimensions see [53]). This 
function is defined on a set of points and is parameterized 
by two values at each such node, in the form Ck + a 2 dk, to 
reflect the a 2 -scaling of our action. The parameters Ck and 
dk are obtained by minimizing the corresponding x 2 . The 
systematic error of the a — > limit is determined by vary- 
ing the node positions. We find that lattice discretization 
errors become large at high magnetic fields due to satu- 
ration of the lattice magnetic flux [44], therefore we only 
include points with Nb/N 2 < 0.1. In Fig. 1 we also show 
the continuum limit of the difference A(£„ + T, d ) jl. 

Next, we address the condensate at nonzero temperature, 
carrying out a similar continuum extrapolation for AS as 
at T = 0, using three lattice spacings with N t = 6, 8 and 
10. The increase of the difference AE(-B) is qualitatively 
similar for zero and nonzero temperatures in ^PT and in 
the PNJL model (see below). In QCD, however, the sit- 
uation is quite different: in Fig. 2 we plot the continuum 
extrapolated lattice results for A(X M + £ d )/2 as functions 
of B for several temperatures, ranging from T = up to 
T = 176 MeV. Note that the transition temperature varies 
from T c (eB = 0) « 158 MeV down to T c (0.9 GeV 2 ) « 
138 MeV [44]. The increasing behavior of AY,(B) at low 
temperatures (T < 130 MeV) corresponding to magnetic 
catalysis continuously transforms into a hump-like struc- 
ture in the crossover region (T = 148 MeV, 153 MeV) 
and then on to a monotonously decreasing dependence 
(T > 163 MeV). We remark that — although in the high 
temperature limit the condensate and its dependence on B 
are suppressed — at T > 190 MeV A£(£?) again starts to 
increase. Furthermore, we note that the strange condensate 
AE S (with a definition similar to that in Eq. (4)) does not 
exhibit this complex dependence on B and T but simply 
increases with growing B for all temperatures. This shows 
that the partly decreasing behavior near the crossover re- 
gion only appears for quark masses below a certain thresh- 



old m thr , inbetween the physical light and strange quark 
masses, m ud < m t hr < m s- 

Comparison to effective theories/models. — In Fig. 3 we 
compare our zero temperature QCD result for A(£„ + 
S d )/2 as a function of B to the %PT prediction [13-15, 54] 
and to that of the PNJL model [18, 55], both at physical 
pion mass. We see that the %PT prediction describes the 
lattice results well up to eB = 0.1 GeV 2 , while the PNJL 
model works quantitatively well up to eB = 0.3 GeV 2 . 
Note that, since the Polyakov loop at zero temperature van- 
ishes, in the limit T — > the PNJL model becomes indis- 
tinguishable from the NJL model with the same couplings. 

In Fig. 4, the condensate Eq. (4) as a function of T is 
compared to %PT and to the PNJL model for different mag- 
netic fields. At B = we use the continuum extrapolation 
for the condensate presented in Ref. [50] (where lattices up 
to N t = 16 were employed), and complement this with the 
differences AT,(B) shown in Fig. 2. In addition to the con- 
tinuum extrapolated lattice data we plot the ^PT curves for 
B = [35] and for B > [14, 15, 54], together with the 
PNJL model predictions [18, 55]. The results indicate that 
XPT is reliable for small temperatures and small magnetic 
fields, eB < 0.1 GeV 2 , T < 100 MeV. (We remark that 
the inclusion of the hadron resonance gas contribution to 
the condensate in x?T [35] improves the agreement with 
lattice results, as was shown at B = in Ref. [50]. One 
would expect a similar improvement at B > 0.) Since 
the PNJL model condensate is calculated using a Polyakov 
loop effective potential that was obtained from Nj = 2 
lattice results [18], differences between the model and our 
Nf = 1 + 1 + 1 results at T > are expected to be large, as 
both the transition temperature and the transition strength 
(the slope of the condensate at T c ) strongly depend on the 
number of flavors. To enable a comparison, we linearly 
rescaled the temperature axis (only for the PNJL curves) 
to match our lattice inflection point at B = 0. Never- 
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theless, the Independence of the condensate for the PNJL 
model also reveals qualitative differences in comparison to 
the QCD results. 

Finally, in Fig. 5 we plot AS U — AT, d = S u — S d 
as a function of the temperature for several magnetic field 
strengths. At zero magnetic field isospin symmetry is exact 
since we employed mass-degenerate light quarks. As B in- 
creases, due to the difference between the electric charges, 
Yi u — T, d develops a temperature-dependence similar to that 
of (S M + T ld )/2, see Fig. 4. The results for T, u ± £ d are 
also listed in Table I. 

Summary. — We determined the QCD light quark con- 
densates at nonzero external magnetic field strengths for 
physical quark masses in the continuum limit. Our results 
are in quantitative agreement with chiral perturbation the- 
ory and PNJL model predictions for small magnetic fields 
and at small temperatures. Note that the constants within 
these parameterizations have not been adjusted to our data 
but were taken from the literature where they have been 
obtained at vanishing magnetic field. Unsurprisingly, yPT 
fails in regions where pions cease to be the essential low en- 
ergy degrees of freedom. While in the hadronic phase low 
energy models qualitatively reproduce the S-dependence 
of the lattice data, they miss an important feature which 
becomes dominant for light quark masses and for tempera- 
tures around T c , see Fig. 2. Clearly, the coupling between 
the magnetic field and the gauge background is enhanced 
near the chiral limit: the smaller the quark mass, the more 
the fluctuations of the gauge field influence the quark de- 
terminant. Thus, for light quarks the indirect interaction 
between the gluonic degrees of freedom and the external 
field becomes more important. A possibility to separate 
this indirect effect would be to consider the sea and va- 
lence contributions to the condensate, as was performed in 
Ref. [39], which we plan to discuss in a forthcoming study. 
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TABLE I. Continuum extrapolated lattice results for the light 
condensates, as functions of T and eB. Columns labeled '+/2' 
contain the average (S u + T, d )/2, while those with ' contain 
the difference E M — E^. Note that the uncertainty of the lattice 
scale gives rise to errors of about 2% in the temperatures. 
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